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ABSTRACT 



, Aims. We recently performed a study of a sample of relaxed, cooling core galaxy clusters with deep Very Large Array observations at 

1.4 GHz. We find that in the central regions of A 1835, A2029, and Ophiuchus the dominant radio galaxy is surrounded by a diffuse 
low-brightness radio emission that takes the form of a mini-halo. Here we present the results of the analysis of the extended diffuse 
radio emission in these mini-halos. 

^ ^, Methods. In order to investigate the morphological properties of the diffuse radio emission in clusters of galaxies we propose to 

I ■ fit their azimuthally averaged brightness profile with an exponential, obtaining the central brightness and the e-folding radius from 

O ' which the radio emissivity can be calculated. We investigate the radio properties of the mini-halos in A1835, A2029, and Ophiuchus 

in comparison with the radio properties of a representative sample of mini-halos and halos already known in the literature. 
C/3 , Results. We find that radio halos can have quite different length-scales but their emissivity is remarkably similar from one halo to the 

other. In contrast, mini-halos span a wide range of radio emissivity. Some of them, like the Perseus mini-halos, are characterized by a 
radio emissivity which is more than 100 times greater than that of radio halos. On the other hand, the new mini-halos in cooling core 
clusters analyzed in this work, namely A2029, Ophiuchus, and A1835, have a radio emissivity which is much more typical of halos 
^ , in merging clusters rather than similar to that of the other mini-halos previously known. 

CO , Key words. Galaxies:clusters:individual: A1835, A2029, Ophiuchus - radio continuum: galaxies 

0^ 

• ■ 1 . Introduction large-scale halos, have a steep spectrum and a very low surface 

' _ brightness. 

O ■ Radio halos are diffuse, low surface brightness, steep-spectrunll| ^^^^ ^^^^^ 30 ^^^^^ j^^j^^ y.^^^^ Giovannini & 

a\ sources (a > 1), permeating the central regions of galaxy clus- pgj.gjjj 2OOO, Bacchi et al. 2003, Govoni et al. 2001, Venturi et 

^ ters. These radio sources are extended on megaparsec scales and 2007, 2008, Giovannini et al. in preparation). Due to their ex- 

;> , are produced by synchrotron radiation of relativistic electrons ^^^^^^^ ^^^.f^^^ brightness and large angular extent (> 10' at 

with energies of- 10 GeVm magnetic fields with B - 0.5-1 //G aredshiftz <0.1) radio halos are best studied at low spatial res- 

(see e.g. Feretti & Giovannini 2008, FeiTari et al. 2008, and ref- Qj^jj^n. In fact, several radio halos were detected by Giovannini 

^ erences therein for recent reviews). (^999^ NRAO VLA Sky Survey (NVSS; Condon et 

Radio halos are typically found in clusters which show al. 1998) and by Kempner & Sarazin (2001) in the Westerbork 

significant evidence for an ongoing merger (e.g. Buote 2001, Northern Sky Survey (WENSS; Rengelink et al. 1997), where 

Govoni et al. 2004). It has been proposed that recent cluster the relatively large beam of these surveys provided the necessary 

mergers may play an important role in the re-acceleration of the sensitivity to large-scale emission to spot these elusive sources, 

radio-emitting relativistic particles, thus providing the energy Dug to the combination of their small angular size and the 

that powers these extended sources (e.g. Brunetti et al. 2001, strong radio emission of the central radio galaxy, the detection 

Petrosian 2001). A major merger event is expected to disrupt of a mini-halo requires a much higher dynamic range and resolu- 

cooling cores and create disturbances readily visible in an X- jion jhan provided by available surveys, thus compUcating their 

ray image of the cluster Therefore, the merger scenario predicts detection. As a consequence, our current observational knowl- 

the absence of large-scale radio halos in symmetric cooling core edge on mini-halos is limited to only a handful of well-studied 

clusters. However, a few cooling core clusters show the presence clusters (e.g., Perseus: Burns et al. 1992; A2390: Bacchi et al. 

of a diffuse synchrotron emission that extends quite far from 2003; RXJ1347.5-1 145: Gitti et al. 2007), and their origin and 

the dominant radio galaxy at the cluster center, forming what physical properties are still poorly understood, 

is called a mini-halo. These diffuse radio sources are extended ^o search for new extended diffuse radio emission in relaxed 

on a moderate scale (typically - 500 kpc) and, in common with ^nd cooling core galaxy clusters, we recendy performed a study 

of a small sample of clusters with deep Very Large Array obser- 

Send offprint requests to: M .Murgia, email m.murgia@ira.inaf.it vations at 1.4 GHz. We find that in the central regions of A1835, 
' We use the convention 5 ^ « v A2029, and Ophiuchus the dominant radio galaxy is surrounded 



2 



M. Murgia et al.: Comparative analysis of the diffuse radio emission in the galaxy clusters A1835, A2029, and Ophiuchus 



Table 1. The mini-halos analyzed in this work and the basic 
properties of the radio images at 1.4 GHz. 



Cluster 


z 


kpc/" 


beam 


(T 










mjy/beam 


A1835 


0.2532 


3.91 


53"x53" 


0.025 


A2029 


0.0765 


1.43 


53"x53" 


0.031 


Ophiuchus 


0.028 


0.55 


92"x92" 


0.16 



Col. 3: Angular to linear size conversion factor; 
Col. 4: FWHM beam of the radio images; 
Col. 5: Image sensitivity 



by a diffuse low-brightness radio emission that takes the form 
of a mini-halo (Govoni et al. 2009, in press; hereafter Paper I). 
Here we present the analysis of these three new mini-halos and 
we investigate their radio properties in comparison with those of 
a statistically significant sample of mini-halos and halos already 
known in the literature and for which good VLA radio images 
at 1.4 GHz are available. In §2 we present an analysis of the 
mini-halos found in A1835, A2029, and Ophiuchus. In particu- 
lar we fit their azimuthally averaged brightness profile with an 
exponential. Through the study of their averaged brightness pro- 
files, we analyze some of their morphological and physical radio 
properties (i.e., length-scale, central brightness, and emissivity). 
In §3 we compare the properties of the mini-halos in A1835, 
A2029, and Ophiuchus with those of other halos and mini-halos 
already available in the literature. In §4 we test the reliability of 
the fitting procedure. The discussion is presented in §5 and our 
conclusions are reported in §6. 

Throughout this paper we assume a ACDM cosmology with 
Hq = 71 km s-'Mpc^^ 0„, = 0.27, and Q.a = 0.73. 

2. Analysis of the mini-halos in A1835, A202g and 
Ophiuchus 

In order to investigate the presence of diffuse extended emis- 
sion in relaxed systems, we recently examined a small sample of 
cooling core clusters with deep Very Large Array observations at 
1 .4 GHz. The combination of resolution and sensitivity of these 
observations allow us to increase the number of known large 
scale diffuse radio emissions in cooling core clusters. Indeed, 
we found that in the central regions of A1835, A2029, and 
Ophiuchus the dominant radio galaxy is surrounded by a dif- 
fuse low-brightness radio emission that takes the form of a mini- 
halo. The basic properties of these clusters are reported in Table 
[1] while the details of the observations are described in Paper I. 
Here we present a morphological analysis of these mini-halos. 
For this purpose, the radio images have been convolved with a 
circular Gaussian beam and corrected for the primary beam at- 
tenuation. 

The brightness of halo-like diffuse radio sources in clusters 
of galaxies decreases with increasing distance from the cluster 
center, eventually falling below the noise level of the radio im- 
ages. The signal-to-noise ratio of the VLA images analyzed in 
Paper I allows us to trace the mini-halo extension out to a dis- 
tance from the cluster center of ~450 kpc, ~250 kpc, and ~200 
kpc for A1835, A2029, and Ophiuchus, respectively (see left- 
panels of Fig. [T]|. However, the size of the diffuse low-surface 
brightness emissions as derived from the 3cr-isophotes seen in 
the radio images may be strongly affected by the sensitivity of 
the radio observations. It is desirable to define a quantity which 
is independent of the signal-to-noise ratio of the radio images 



in order to obtain an unbiased estimate of the size of mini-halos 
and halos. 

Although deviations of the diffuse emission from spherical 
symmetry are often observed, the azimuthally averaged radial 
profiles are indeed quite stable. For this reason, following the 
same approach used by Orrii et al. (2007), we fit the azimuthally 
averaged brightness profiles with an exponential law of the form: 

I(r) = ke-''^' (1) 

where the two independent parameters Iq and are the central 
radio surface brightness and the e-folding radius i.e., the radius 
at which the brightness drops to /o/e. As we will show, the quan- 
tity re can be used to define a length-scale which is relatively in- 
dependent on the sensitivity of the radio images. The exponential 
model is attractive in its simplicity and involves a minimal set of 
free parameters. 

We note that in mini-halos the presence of a strong central 
radio galaxy, which appears as a point source at the low reso- 
lution of our observations, complicates the analysis. Therefore 
its emission must be known accurately in order to separate its 
contribution from that of the mini-halo. 

In Fig. [U we present the resulting azimuthally averaged 
brightness profiles of A1835, A2029 and Ophiuchus, respec- 
tively. The profiles have been obtained by averaging the radio 
brightness in concentric annuli as shown in the left panels of Fig. 
[1] The annuli are as wide as the half FWHM beam. In the statis- 
tical analysis all the discrete sources (except the central source) 
have been excluded by masking them out. Moreover, we con- 
sidered only data points above the 3-sigma radio noise level as 
indicated by the horizontal dashed-dotted line in the plots. 

In order to carefully separate the contribution of the mini- 
halo from that of the central radio galaxy, we fitted the total 
brightness profiles with a central point source plus the radio 
mini-halo diffuse emission 

I{r) ^ Ipsir) + iMH(r). (2) 

In particular the profile of the central point source has been char- 
acterized by a Gaussian of the form 

/ps(r)=/o„^^-<^'^'"™' (3) 

while the brightness profile of the mini-halo has been character- 
ized by an exponential law of the form 

/MH(r) = /„„,e-'■/'■^ (4) 

The resolution of the observations, the masked regions, and 
the sampling of the radial profile in annuli, are all factors that can 
affect the estimate of the best fit parameters and thus we decided 
to include their effects directly in the modelling. 

In order to properly take into account the resolution, the ex- 
ponential model in Eq.|4]is first calculated in a 2-dimensional 
image with the same pixel size and field of view as that ob- 
served, and then convolved with the same beam by means of a 
Fast Fourier Transform. The point-source model in Eq|3]is added 
and the resulting image is masked exactly in the same regions as 
for the observations. Finally, the model is azimuthally averaged 
with the same set of annuli used to obtain the observed radial 
profile. All these functions are performed at each step during the 
fit procedure. As a result, the values of the central brightness, 
Iq, and e-folding radius, provided by the fit are deconvolved 
quantities and their estimate includes all the uncertainties related 
to the masked regions and to the sampling of the radial profile in 
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Fig. 1. The azimuthally averaged brightness profiles of the radio emission in A1835 (top), A2029 (middle), and Ophiuchus (bottom). 
The profiles have been calculated in concentric annuli, as shown in the left panels. All the discrete sources (except the central source) 
have been masked. For each cluster the horizontal dashed-dotted line indicates the 3cr noise level of the radio image. In our analysis 
we considered all the data points above the 3cr noise level. The black line indicates the best fit profile described by an exponential 
law (dashed line) representing the mini-halo emission, and by a central Gaussian profile (dotted line) representing the central point 
source. 
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Table 2. Radio properties of the mini halos derived from the fit procedure 



Cluster 


FWHM 

arcsec 


re 

arcsec 


kpc 


yujy/arcsec^ 


S MH 

mJy 


//Jy/arcsec- 


FWHMps 
arcsec 


S PS 
mJy 


erg s^'Hz^'cm"^^ 


Xred 


A1835 


53 




102^™ 




"■"-0.7 


8.68:0^] 


57+2 


9+0.8 
oz,.z,„7 


3.3!j;2 X 10-« 


0.03 


A2029 


53 


37!* 


531 


9 7^+0.92 
^•'-^-0.64 


18.8:;;3 


129!^ 


57+1 


480.0!|« 


5.4!|;^ X 10-« 


0.5 


OPHIUCHUS 


92 


191+23 


105!1] 


0.58!«:«? 


106.4+^0/ 


7 47+0.22 


102+^ 


29.0+2-;! 


4.7!«;9 X 10-« 


2.9 



Col. 2: FWHM of the circular Gaussian beam; 

Col. 6: Mini-halo radio flux density calculated up to Sr^; 

Col. 10: Average radio emissivity over the volume of a sphere of radius 3;'^., /:-coiTected with = 1. 



annuli of finite width. The fit procedure has been implemented 
in the software FARADAY (Murgia et al. 2004). 

The best fit model is shown as a continuous black line in the 
right panels of Fig. [T] The mini-halo contribution is indicated 
by the green dashed line, while the contribution of the central 
point source is indicated with the red dotted line. The best fit 
values for both the mini-halos and the central components are 
reported in Table |2] together with their respectively errors and 
the coiTesponding;^^^. 

Overall, the observed profiles are well fit by the exponen- 
tial model. In A2029 and Ophiuchus the emission of the mini- 
halo extends for several beams from the central point source. 
It is important to note the fundamental difference between the 
method based on the measurement of rmax, the mini-halo max- 
imum size from the 3cr-isophote above the nois^E and the fit 
of the e-folding radius, r^. In fact, the mini-halo in A2029 has 
a larger extension with respect to that in Ophiuchus: the two 
mini-halos have r„,ax ~ 250 kpc and r,„ax ~ 200 kpc, respec- 
tively. However, the mini-halo emission in A2029 falls off more 
rapidly and the corresponding e-folding radius = 53!g kpc is 
much smaller than that in Ophiuchus, = 105^] j kpc. The fit is 
critical for A1835 (see §4). This is the most distant cluster of the 
three and the emission of the mini-halo is heavily blended with 
that of the central point source. The maximum extent for A1835 
is as large as r^ax ~ 450 kpc, while the e-folding radius obtained 
by the model fit results, with a large uncertainty, is = 105^2" 
kpc. 

The central brightness ranges from /q = 2.75 /iJy/arcsec^ 
for A2029 to Iq = 0.58 /iJy/arcsec^ for Ophiuchus, indicating 
a slightly tendency for the smaller mini-halos to have a higher 
central brightness (see §3). In all three cases the best fit model 
results in a central Gaussian with FWHMps - crp^ x 2.35 > 
FWHMbeanu as cxpcctcd for a point like or a slightly resolved 
source. From Iq and we then calculated the flux density S mh 
(see Table |2]i of the mini-halos. In particular, the flux density is 
derived by integrating the brightness profile I{r) up to a radius r' 
from the cluster center: 

Smh ^2n { I(r)rdr = 27r[l + e~'-'''-'(-r' /r, - 1)] ■ r%. (5) 
Jo 

We calculated the flux density of the mini-halos by integrat- 
ing the surface brightness of the best fit exponential model up to 
r' = 3re, giving: 

SMH^27Tf-r% (Jy) (6) 

2 In this paper, we define r,„n.v as the distance from the cluster cen- 
ter at which the azimuthally averaged brightness profile of the diffuse 
emission reaches the 3cr noise level of the radio image. 



where and /q are in units of arcsec and Jy/arcsec^, respec- 
tively, while / = (1 - ^ 0.8. Thus, the flux density in Eq|6] 
consists of about 80% of the total model flux that would be ob- 
tained by extrapolating r' to +oo. It must be noted that the model 
flux calculated from Eq|6]may differ from the flux density mea- 
sured by integrating the radio brightness up to the 3cr-isophote. 
For instance, in the case of Ophiuchus r^ax < 3rj, and the model 
flux of 106 mJy is higher than the flux density of about 80 mJy 
which is obtained by integrating up to the 3cr-isophote. Instead, 
in the case of A2029, r^ax > Sr^ and the model flux density of 
18.8 mJy is lower than the flux density of 22.5 mJy obtained by 
integrating up to the 3cr-isophote. 

In column 9 of Tab. 2, we also provide the flux density of the 
central point source derived from the fit. The flux density of the 
central point source, 5 ps , is calculated through: 

^''^'8^)'^^"^''^''' ^^^^ 

where FWHMps is measured in arcsec and /q^j in Jy/arcsec^. 
Since mini-halos are optically thin sources, the radio brightness 
/, at a projected distance r from the cluster centre, is related to 
the radio emissivity J{R) through the line-of-sight integral: 

Kr)^^ f J(R)dl (8) 

J source 

where R - Vr^ -H P. 

We calculated a volume-averaged radio emissivity by sup- 
posing that all the flux density in Eq.|6] comes from a sphere of 
radius 3re: 

{J} ^7.7xl0-'^\l+zf^" ■- (ergs^'cm^^Hz"') (9) 

re 

where and /q are in units of kpc and //Jy/arcsec^, respectively. 
The factor (1 -H z)^'^" takes into account of both the A;-correction 
and cosmological dimming of the surface brightness with red- 
shift. The radio emissivity for the three mini-halos (A:-coiTected 
with a spectral index a = 1) is reported in column 10 of Tab. 2. 

3. Mini-halos and halos comparison 

Merging and cooling core clusters containing extended diffuse 
radio sources (halos and mini-halos respectively) show different 
X-ray properties, indicating a different evolutionary state. A sys- 
tematic and homogeneous comparison between the properties of 
extended diffuse radio sources in merging systems with those 
in cooling core clusters, may help in understanding if the energy 
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necessary to trigger their radio emission may be related to differ- 
ent physical processes. Therefore, we analyzed the azimuthally 
averaged brightness profiles in a sample of 12 clusters contain- 
ing a central radio halo, previously imaged by us. Moreover, in 
order to investigate also a significant number of mini-halos, we 
reduced and analyzed the VLA data at 1.4 GHz of the mini- 
halos in RXJ1347.5-1145 (Gitti et al. 2007), A2390 (Bacchi et 
al. 2003), and Perseus (Pedlar et al. 1990), previously reported 
in the literature. The cluster list, together with their references is 
given in Table [3] 

Indeed the analysis lacks some well known radio halos and 
mini-halos, however, the number of clusters analyzed here is rep- 
resentative of these classes of radio sources and contains both ex- 
tended (e.g. A2163, A2744) and small halos (e.g. A401, A2218) 
in merging clusters and several mini-halos in cooling core clus- 
ters. In Fig. |2] and [3] we present the resulting azimuthally aver- 
aged brightness profiles of the radio halos, while in Fig. |4] we 
present those of the three mini-halos. The profiles have been ob- 
tained by averaging the radio brightness in concentric annuli as 
shown in the inset of each panel. 

The analysis of the three previously known mini-halos has 
been performed with the same approach proposed for A1835, 
A2029, and Ophiuchus - by fitting the brightness profiles with a 
central Gaussian plus an exponential law, and by masking all the 
discrete sources, except the central one. The analysis of the halos 
has been performed by fitting the brightness profiles with only 
an exponential law and by masking all the discrete sources. As 
in the case of A1835, A2029, and Ophiuchus the diffuse syn- 
chrotron emission both of the halos and mini-halos are quite 
well fit by an exponential law characterized by a central surface 
brightness /q and an e-folding radius r^. The best fit values of 
the fit procedure, together with their respectively errors, and the 
corresponding are reported in Tab. 3. 

In the top panel of Fig. |5]we show the best fit central bright- 
ness /o (yuJy/arcsec^) versus the length-scale (kpc) both for ha- 
los (blue triangles) and mini-halos (red dots). The length scale r^, 
in this sample of halos ranges between ^70-400 kpc. While the 
central brightness Iq ranges between ^0.4 - 3 //Jy/arcsec^. All the 
halo sources of the sample populate a well defined region of the 
Iq — re plane. There is a tendency for halo sources with a higher 
central brightness to have a greater length-scale. By comparing 
the location of the halo sources in the Iq - plane with the dot- 
ted lines of constant emissivity we found, for the first time, that 
all halos have a markedly similar radio emissivity 0. 

In contrast, the mini-halo distribution in the Iq - plane 
appears much more scattered. In our sample there are no mini- 
halos with a Ve larger than ^ 100 kpc. They have a tendency 
to populate the low Vg, while the central emissivity can be two 
orders of magnitude higher than in halos. By comparing the lo- 
cation of the mini-halo sources in the Iq - plane with the lines 
of constant emissivity we found that although some mini-halos 
(like Ophiuchus and A1835) are quite comparable to the halo 
sources, in general mini-halos appear clearly separated from the 
halo sources, showing a larger spread in radio emissivity. 

As previously noted, there is a tendency for the smaller mini- 
halos of our sample to have a higher central brightness. But this 
trend must be considered carefully, because a selection effect 
may be present. The dashed line in top panel of Fig. |5]indicates 



^ However, we cannot exclude the possibility that large halos with 
faint surface brightness (/o < 0.5 pjy/arcsec-) could have been missed 
by the current searches at 1.4 GHz such that we are just seeing the "tip 
of the iceberg". As an example, the 3cr sensitivity level of the NVSS is 
indicated by the dot-dashed line in the top panel of Fig. 5. 



the minimum central brightness required for a mini-halo (or a 
halo) with a given e-folding radius to be imaged with at least 4 
resolution elements (beams) in a deep image with a sensitivity of 
cr -25 yuJy/beam, and a 25 arcsec beam. These numbers are ap- 
propriate for the deepest images taken with the VLA at this res- 
olution and frequency (see e.g. Govoni et al. 2005). The detec- 
tion limit is calculated assuming a putative redshift of z = 0.18 
(the average redshift of the mini-halos and halos analyzed in this 
work). In practice, for a mini-halo (or a halo) to be detected, it 
is required that r,„a^ > 2 ■ FWHMteam, where r„u,x is implicitly 
defined by the equation: lQexp(-rmax/re) - 3cr which implies 
/() > 3crexp(2 ■ FWHMheam/fe). The shaded region limited by 
this detection threshold is hardly accessible by current radio in- 
terferometers at 1.4 GHz. For instance, an object with rg ^ 50 
kpc would be detectable only if brighter than Iq > l/iJy/arcsec^. 
As a result, the observed tendency for the smaller mini-halos to 
be brighter must be considered prudently. In any case, the possi- 
ble selection effect mentioned above seems not to influence the 
result that mini-halos have systematically lower values of than 
halos. Mini-halos appear effectively smaller than halos and this 
could reflect the higher density of the intergalactic gas in the 
central regions of cooling clusters. 

For a better comparison between halos and mini-halos in 
terms of emissivity, in Fig. |5] (bottom panel) we show the dis- 
tribution of the radio emissivity (calculated by using Eq|9]l both 
for halos (in blue) and mini-halos (in red). We find that radio 
halos have a surprisingly narrow radio emissivity distribution 
around an average value of about lO"'*^ erg s^'cm^^Hz"'. In con- 
trast, mini-halos span a wide range of radio emissivity. Some of 
them, like the Perseus mini-halo, would be characterized by a 
radio emissivity more than 100 times larger than that of a typical 
halo, quite in agreement with the recent finding by Cassano et 
al. (2008). We note, however, that the three new mini-halos an- 
alyzed in this work are all characterized by brightness and size 
very similar to that of the smaller halos and thus their radio emis- 
sivity is comparable to that of halos in general. 

As a further check, we also examined the trend of the ra- 
dio power (calculated from the flux density in Eq|6]l versus the 
e-folding radius of mini-halo and halos, see Fig. |6] . This plot 
is similar to that presented in Cassano et al. (2008) except that 
they calculated the radius as Vm/n ^ ''max, rmm and r„,a.i being the 
minimum and the maximum radius measured using the 3cr radio 
isophotes. As pointed out by Cassano et al. (2007), the power- 
ful radio halos are the more extended sources while the fainter 
ones are smaller in size. Fig.|6]shows indeed a tight coiTelation 
between the radio power and the e-folding radius of radio ha- 
los. We find that the radio power increases as Pi 4 oc r^. In our 
analysis, the slope of the correlation is in agreement with the 
finding that the radio emissivity has a very small scatter among 
radio halos. In fact, the radio power scales as f 1 4 - J V, where 
V is the source volume and J is the radio emissivity. If we as- 
sume that the emissivity does not change from one halo to the 
other it follows that Pi 4 oc r^. The dotted reference line in Fig.|6] 
represents the expected radio power if we assume for the emis- 
sivity the average value we found for radio halos {J) ^ 10""^^ 
erg s 'cm^^Hz ' and for the volume V the volume of a sphere a 
radius 3r,.: log P1.4 ^ 23.52 + 3 log(rJlOOkpc). 

We note that Cassano et al. (2008) found also a coiTela- 
tion between the power and the size of mini-halos. In fact, the 
three mini-halos we have in common (A2390, RXCJ1347, and 
Perseus) all seem to be aligned in a similar but offset correlation 
with respect to radio halos. However, with the addition of the 
new objects found in Paper I and analyzed here, we cannot de- 
fine any clear coiTelation for the mini-halos, in agreement with 
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Fig. 2. Azimuthally averaged brightness profiles of the radio halos in A2744 (Govoni et al. 2001), A665 (Giovannini & Feretti 
2000), A2219 (Bacchi et al. 2003), A2255 (Govoni et al. 2005), A773 (Govoni et al. 2001), A545 (Bacchi et al. 2003). The profiles 
have been calculated in concentric annuli, as shown in the inset panels. All the discrete sources have been masked. For each cluster 
the horizontal dashed-dotted line indicates the 3cr noise level of the radio image. In the analysis we considered all the data points 
above the 3cr radio noise level. The black line indicates the best fit exponential profile. 
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Fig. 3. Azimuthally averaged brightness profiles of the radio halos in A2319 (Feretti et al. 1997), A2218 (Giovannini & Feretti 
2000), A2163 (Feretti et al. 2001), A401 (Bacchi et al. 2003), A2254 (Govoni et al. 2001), RXCJ1314.4-2515 (Feretti et al. 2005). 
The profiles have been calculated in concentric annuli, as shown in the inset panels. All the discrete sources have been masked. For 
each cluster the horizontal dashed-dotted line indicates the 3cr noise level of the radio image. In the analysis we considered all the 
data points above the 3cr radio noise level. The black line indicates the best fit exponential profile. 
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the literature (Gitti et al. 2007, Bacchi et al. 2003, Pedlar et al. 1990). The profiles have been calculated in concentric annuli, as 
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level. The black line indicates the best fit profile described by an exponential law (dashed line), representing the mini-halo emission, 
and by a central Gaussian profile (dotted line), representing the central point source. 



the large spread in emissivity shown in Fig. |5] The mini-halos 
in A1835, Ophiuchus, and A2029 appear clearly separated from 
the other mini-halos although this separation into two distinct 
groups is likely just due to the small-number statistics. 

As a final remark, we stress once again that our analysis 
should be relatively independent of the noise level of the radio 
images. To ensure that our approach is not dependent on the sen- 
sitivity of the radio observation we retrieved for A2744 (one of 
the clusters analyzed here) the NVSS image. We derived the az- 
imuthally averaged brightness profile as computed for the deeper 
observations (same rings and mask), and we then fit the NVSS 
profile with an exponential law. The resulting fit is shown in the 
left panel of Fig.|2l The contours of the NVSS image are shown 
on the top right panel of Fig.|7]in comparison with the contours 
of the deeper VLA images on the bottom right panel. As ex- 
pected, the radius as derived from the contour levels yields very 
different values in images with different sensitivities. One should 
consider that the same effect may appear in images with compa- 
rable sensitivity of objects characterized by the same intrinsic 



size but different brightness. Fainter halos (or mini-halos) tend 
to "sink" into the noise appearing smaller than they really are. 
On the other hand, in comparing the fitting results of the NVSS 
profiles in Fig. |7] with the fitting results of the deeper images in 
Tab. [3] we find the results to be in good agreement. Therefore 
we are confident that the method discussed here, based on the 
determination of length-scale and central brightness in extended 
structures, is less subject to the S/N ratio of the radio images and 
thus better suited for comparative studies of the general proper- 
ties of diffuse halo-like sources. 



4. Monte Carlo statistical analysis of the fitting 
procedure 

The fit of the exponential model described in §2 could be crit- 
ical in the case of those mini-halos in which the e-folding ra- 
dius is smaller than the FWHM of the central point source. 
The mini-halos in A1835, A2029, RXCJ1347, and A2390 are 
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all in this category. However, it should be noted that not only 
the observing resolution but also the brightness of the mini-halo 
is crucial for the determination of the e-folding radius. If the 
mini-halo is bright enough, the exponential profile will emerge 
from the central point source and it will be observable out to 
a large distance from the cluster center, making the fit of the 
e-folding radius r^, possible. The example of A2029 is represen- 
tative of this situation. The e-folding radius derived from the fit, 
re = 37" + 4", is below the FWHM beam of the observation 
which is 53". Nevertheless, the relative intensity of the mini- 
halo with respect to the noise permits the fit of the exponential 
model from the external part of the mini-halo profile that does 
not overlap with the central point source. The limiting cases of 
A1835 and RXCJ1347 (the two more distant mini-halos in our 
sample) are more critical and the behaviour of the fitting proce- 
dure must be carefully checked. 

In order to assess the reliability of the best fit parameters, 
we created a set of Monte Carlo simulations with the aid of the 
Cybersar-OAC computer cluster. For each mini-halo we varied 
the values of /q and r^, in a regular grid that contained the best fit 
parameters reported in Tab. |2] and Tab. [3] For each pair of Iq and 
re values, we created a synthetic image with the same field of 
view, pixel size, FWHM beam, noise level, and masked regions 
as that observed. A central point source with the same intensity 
as the observed one is added to the exponential model. An ex- 
ample of such image, drawn from the Monte Carlo simulation of 
the Ophiuchus mini-halo, is presented in left panel of Fig.|8] We 
then feed the synthetic image into the fit procedure and check 
for any systematic difference between the input and the output 
values of Iq and rg. In the right panel of Fig. [8] we show the fit 
obtained from the synthetic image on the left. We can see that 
in the case of the Ophiuchus mini-halo, which is particularly ex- 
tended, the difference between input and the output values of Iq 



and is fully consistent with the uncertainties reported in Table 
|2] In Fig.|9]we present the result of the Monte Carlo simulations 
for all six mini-halos analyzed in this work. For each mini-halo 
we plot the grid of the input values of /q and r^ (crosses) with the 
relative arrows representing the displacement of the output pa- 
rameters from the fit procedure. The dashed line represents the 
detection limit for each mini-halo and has been traced with the 
same criterion as in top panel of Fig. |5] Here the exact values 
for the FWHM beam and noise level of the specific observations 
have been used. In each panel, the dot indicates the values of the 
actual best fit parameters with their uncertainties. Clearly, for 
bright and extended mini-halos, the fit procedure is able to re- 
cover the original parameters with high precisions. Obviously, 
near the detection limit the situation can be much worse and 
large systematic effects are observed. However, we note that the 
fit procedure is able to recover the input e-folding even if it is 
significantly smaller then the FWHM beam, provided that the 
value of Iq is sufficiently high (see e.g. the case of A2390 with 
re - 10" and Iq > 10 //Jy/beam). 

In summary, the results of the Monte Carlo simulation tell 
us that the fit procedure is able to constraint the model param- 
eters with accuracy and without systematic effects in the cases 
of Ophiuchus, Perseus, A2029, and A2390. The mini-halos in 
A1835 and RXCJ1347 are very close to the detection limit and 
the best fit parameters could be affected by a systematic bias. 
In particular, they could be smaller and brighter than what is 
inferred by the fit to the exponential model, although the bias 
appears to be comparable to the reported uncertainties. 

5. Discussion 

The exponential model fitting method provides an estimate of 
the brightness and size of the diffuse, halo-like, radio emission 
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Fig. 7. Left: fit of the azimuthally averaged brightness profiles calculated from the NVSS image of the cluster A2744. Right: 
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while the sensitivity of the deeper image is 0. 1 mJy/beam. The radius as derived from the contour levels depends on the sensitivity 
of the radio image, while the fit procedure give similar results in images at diff'erent sensitivity. 



in cluster of galaxies using a minimal set of free parameters. 
In this section we discuss the exponential profile in comparison 
with alternative models which may have a much clearer physical 
meaning, but require the introduction of a number of additional 
assumptions and/or free parameters. 

One main difficulty in the study of the non-thermal radio ha- 
los and mini-halos in clusters of galaxies is that the synchrotron 
emissivity essentially traces the product of electron and mag- 
netic field energy densities so disentangling the two contribu- 
tions is not easy. Another source of uncertainty is related to the 
effective shape of the energy spectrum of the relativistic elec- 
trons across the cluster which is poorly known. Given their com- 
paratively short radiative life-times, the synchrotron electrons 
cannot diffuse into the large volume of space involved but rather 
they must be either continuously injected and/or re-accelerated 
in situ throughout the intra-cluster medium (e.g., Dennison 1980, 
Jaffe 1977, Dolag & Ensshn 2000; Brunetti et. al 2001). The 
spectral energy distributions for the cosmic ray electrons pre- 
dicted by these models can be quite different, however from the 
current data we are not able to discriminate between the two 
processes definitively. Our intent here, is to highlight the role of 
the magnetic field in determining the radial profile of the radio 
emission. For illustrative purposes, in the following we discuss 
simplified models for the distribution of particles and magnetic 
fields in clusters that are often discussed in literature and we 
compare them with the exponential profile. 



5.1. Injection models 

We can compare a continuous injection of cosmic ray electrons 
in the intra-cluster medium with a power law distribution in the 
energy range from y,„i„mec'^ to jmaxnieC^ {jmax » Jmin)- In the 
stationary case, and for an emission spectrum with index a - \, 



from standard synchrotron formulas (Rybicki & Lightman 1979) 
we obtain for the radio emissivity: 

Jy = 2.21 X 10-'° fi^cVG^, (erg s-'cm-^^Hz-')(10) 

"nC "CMB 

where Bcmb - 3.25(1 + z)^ fiG is the equivalent magnetic field 
associated with the cosmic microwave background while Q^i 
is the production rate of relativistic electrons in units of erg 
s 'cm"^. In Eq.[TO]we supposed that the electron population is 
isotropic and we average over all the possible directions between 
the magnetic field and the line-of-sight, i.e. the field is com- 
pletely tangled. 

The energy density of the cosmic ray electrons can be calcu- 
lated as 

- 7.74 X lO^O— ^^y;,|„ (erg cm-^). (11) 

"tiG "CMB 

In the case of a perfect equipartition between the energy den- 
sities of particles and fields (m^,/ - B^/Sn), we have that the radio 
source is near minimum energy conditions and the radio emis- 
sivity at 1.4 GHz is given by: 

7i.4 - 8.12 X m-''y„MQoBlc (erg s-'cm-^Hz"'). (12) 

The synchrotron emissivity is indeed very close to the observed 
value of {J} ^ 10"^^ erg s^'cm^^Hz"' in the case of radio halos 
for a low energy cut-off of 'X,„in = 100 and an average magnetic 
field strength of about 1 juG (changing y„„„ from 10 to 1000 
causes B to vary from 1.8 to 0.5 juG). Assuming a redshift of 
z=0.18 (the average distance of the radio halos in our sample), 
from Eq.[TT] we find that the required energy supply to explain 
the observed radio halos emissivity in a minimum energy con- 
dition is Qei - 1.1 X 10"-" ergs" 'cm"-'. Note that this power is 
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Fig. 8. Synthetic image (left panel) and fit profile (right panel) drawn from the Monte Carlo simulation of the Ophiuchus mini-halo. 



about an order of magnitude larger than the total synchrotron lu- 
minosity and one order of magnitude less than the thermal X-ray 
luminosity. 

The equipartition magnetic field strength quoted above 
should be regarded as volume averages. However, there are ob- 
servational and theoretical arguments to support the idea that in- 
side a cluster the magnetic field strength scales as a function of 
the thermal gas density, n,/,, so: 



B(R) = Bo 



«r/,(0) 



(13) 



where the index rj is expected to be in the range 0.3 - 1 (e.g., 
Dolag et al. 2001, Dolag et al. 2002, Murgia et al. 2004, Guidetti 
et al. 2008, Battaglia et al. 2009). 

If the density profile of the thermal gas is described by 
the standard y6-model (Cavaliere & Fusco-Femiano 1976), from 
Eq.[T2]we have 



(14) 



where rc is the cluster's core radius. The radio brightness 
profile is: 



hAr) = /o(l + rVr.^)-"^"^"-' 



(15) 



where r is the projected distance from the cluster center. It is 
worthwhile to note that the ratio between the radio and X-ray 
brightness is proportional to 



2 , 2x-3/S(2;,-l) 



(16) 



where T is the gas temperature. Thus, if r] - 0.5, a linear 
relation is expected between the radio and the X-ray brightness 
for a perfectly isothermal j6-model. 

We now reconsider the radio halos in A2744 and A665, for 
which an average spectral index a - 1 has been determined 



from VLA radio images at 327 MHz and 1.4 GHz by OiTii 
et al. (2007) and Feretti et al. (2004), respectively. In the top 
panels of Fig. 10, we trace the equipartition model profiles by 
converting to our cosmology the y6-model parameters reported 
by Kempner & David (2004) and Birkinshaw et al. (1991) for 
A2744 (r, = 382 kpc and /3 = 0.79) and A665 (r, = 291 kpc and 
/3 - 0.66), respectively. The equipartition model is very sensi- 
tive to the radial profile of the magnetic field strength. In order to 
show this behaviour, we fix Iq in Eq{T5]and choose three particu- 
lar values for t]. A constant field, rj -Q, results in a profile which 
is too flat with respect to the data while a magnetic field whose 
strength decreases with the distance as the gas density, rj - I, 
results in a profile which is too steep. A very good agreement 
with the data is found for ^ 0.4 - 0.5 for both radio halos, im- 
plying a magnetic fields whose energy density is proportional to 
that of the thermal gas. Indeed, in A2744 a linear correlation be- 
tween the radio and X-ray surface brightness has been observed 
(Govoni et al. 2001). Moreover, a comparison of the equiparti- 
tion model with rj - 0.5 and the exponential fit shows that the 
two profiles are strikingly similar within one megaparsec from 
the cluster center. The differences between the two profiles be- 
come appreciable only at much larger distances, where the radio 
halo emission has fallen below the noise level. 

In the bottom panels of Fig. 10 we show the exponential and 
equipartition profiles for the mini-halos in the Ophiuchus and 
Perseus clusters. We adopt a core radius of rc - 199 kpc and 
P - 0.747 for Ophiuchus, and - 44 kpc and (3 - 0.54 for 
Perseus (Chen et al 2007). We note that the equipartition pro- 
file in Eq[T5] should be considered only outside the cool core 
where the j6-model still provides a reasonable description of the 
gas density. Moreover, Ophiuchus lacks any spectral index in- 
formation while the mini-halo in Perseus has an average spec- 
tral index which is slightly steeper than what we have assumed 
here {a ^ 1.2; Gitti et al. 2004). Apart from these caveats, how- 
ever, the exponential and the equipartition profiles are remark- 
ably similar in both mini-halos. As found for the radio halos 
in A2744 and A665, the mini-halo in Ophiuchus is well de- 
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Fig. 9. Monte Carlo simulations for the six mini-halos analyzed in this work, see text for further details. 
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scribed by an equipartition model with 77 ^ 0.5. In the case of 
the Perseus mini-halo, however, we found 77 ^ 0.8, which im- 
plies a particularly steep magnetic field radial profile (this re- 
sult does not change significantly by assuming a - \.2 instead 
of a = 1). Indeed, Burns et al. (1992) suggested that the ab- 
sence of a cluster-wide magnetic field inhibits a large scale halo 
in Perseus. 

We conclude that the exponential profile is very close to 
the expectations of a simple equipartition model if the magnetic 
field energy density roughly scales as the thermal gas density. In 
the framework of the equipartition model, the physical meaning 
of the 3re length-scale obtained by the exponential fit is that it 
marks the point at which the magnetic field strength is decreased 
to about half the value at the cluster center (see bottom panel of 
Fig. 10). Hence, mini-halos would appear smaller because the 
magnetic field falls off more rapidly with radius. 

It is worthwhile to note that if Z? < Bcmb and 77 ^ 0.5 (as 
for A2744, A665 and Ophiuchus), the equipartition condition 
impUes that the injection rate should scale roughly as the ther- 
mal gas density Q^i oc n,h. In models that consider a continuous 
production of secondary electrons by hadronic collisions in the 
intra-cluster medium (Dennison 1980), the injection rate of cos- 
mic ray electrons is proportional to the product of the densities 
of the gas and cosmic ray protons, Q^i oc n,kncRp- Usually, it is 
assumed that ricRp n,k and as a results Q^i oc n^^^. As a con- 
sequence secondary models predict brightness profiles that are 
generally too steep to explain observed trends in radio halos. 
However, the radio brightness profiles in the Perseus mini-halo 
is steeper and it can be described by secondary models as well 
(Pfrommer & Ensslin 2004). 

5.2. Re-acceleration models 

The observed radio profiles may be also described in the frame- 
work of the re-acceleration scenario (e.g. Gitti et al. 2004). 
An over-simplified re-acceleration model may be considered in 
which the relativistic electrons energy spectrum is essentially 
peaked at a characteristic energy, jmaxmc^, at which the re- 
acceleration gains balance the radiative losses. In the stationary 
case, the energy of the cosmic ray electrons can be calculated as: 

Uei - 7.74 X 10'" 2 ^'L TmL (erg cm-^) (17) 

"nC "CMB 

where Qei = Uei is the energy supply to the re-acceleration pro- 
cess in units of erg s^'cm^^*. 

The emission spectrum will extend up to a maximum fre- 
quency of v,„fl.i. = 4.2 X 10^'* B ±pc7max GHz, beyond which it 
cuts off exponentially. The synchrotron emissivity in c.g.s. units 
is: 

Jy = 2.74 X 10'"' 2 Bl^KLF(y/y,na.) (18) 

"fiG "CMB 

where F(v/vmax) is the usual synchrotron kernel (Pacholczky 
1970) while v and v„ax are expressed in GHz. 

In order to reproduce a spectral index = 1 between 

325 MHz and 1.4 GHz, it should be v„,„ ^ 0.56 ■ ('l -Hz) GHz. At 
a frequency of v = 1.4 GHz and at z = 0.18, the average redshift 
of the sources in our sample, we can calculate F(v/Vmax) - 0.27 
(a steeper spectrum with - 2 would imply v„,ax - 0.3 1 ■ 

(l+z) GHz and F(v/v,„„x) ^ O.dS) . 

In the case of a perfect equipartition between the energy den- 
sities of particles and fields, Ud = B~/En, we find that the radio 



source is near the minimum energy condition and the radio emis- 
sivity at 1.4 GHz is: 

7i.4 ^ 7.65 X 10"'*'B;J^ (ergs-'cm^^Hz-'). (19) 

The observed synchrotron emissivity in the case of radio ha- 
los, (7) ^ lO"'*^ erg s^'cm^-'Hz"', implies an average magnetic 
field strength of about 0.3 /kG. At a redshift z = 0.18, fromEq.fTSl 
it follows that the required energy supply is Q^i - 2.2 x 10"^" erg 
s 'cm"^, about one order of magnitude larger than the required 
power found in the previous section in the case of the continuous 
injection model. 

A general result from Eq.[T7] is that there exists a critical 
value for the magnetic field strength that minimize the accel- 
eration efficiency = Qei/uei needed to radiate at a given v„ax, 
described hy B - Bcmb I V3. In the particular case of the aver- 
age redshift of halos in our sample, the field that minimizes the 
re-acceleration efficiency results in Z? ^ 2.6 /iG. Given this value 
for the magnetic field, from Eq.[T8] it follows that the required 
energy supply to the re-acceleration process is 2e/ - 3.6 x 10"^^ 
erg s^'cm"^, i.e. significantly less than in the equipartition case 
above. However, the energy density of cosmic ray electrons, 
Uei ^ 1.3 X 10"'^ ergcm"^, is much lower than that of the field, 
ub - 2.7 X lO"'-' erg cm"-', and thus the radio source cannot be 
in a condition of minimum energy. 

In Fig. 11 we show the comparison of the re-acceleration 
model versus the exponential fit for the radio halo in A2744. The 
average spectral index of this radio halo is fairly constant with 
a\2SMH~ ~ ^ (Orrii et al. 2007) which justifies the use of the same 
v,„ax at different distances from the cluster centre. By assuming 
that B oc from Eq.[T9]it follows that the radio brightness 
profiles is: 

/l.4(r) = /o(l + ,2/,2)-5-25/S.+0.5^ (20) 

If the equipartition energy condition u^i - ub is met at any 
radius, the re-acceleration model yields a central magnetic field 
of Bo - 0.6 yuG with a scaling index of rj - 0.47. If we fix 
Bo = 3.2 yuG, we minimize the re-acceleration efficiency at the 
cluster centre but the radio halo is out of equipartition since the 
total energy density is dominated by the magnetic field energy 
density (m^,/ - 3 x IO^^mb). Whatever the case, we find that the 
exponential profile is also very close to the expectations of the 
simple re-acceleration model considered here. 

6. Conclusions 

Mini-halos in clusters are still poorly understood sources. They 
are a rare phenomenon, having been found so far only in few 
clusters. A larger number of mini-halos and better information 
on their physical properties will be necessary in order to dis- 
criminate between the different mechanisms suggested for trans- 
ferring energy to the relativistic electrons that power the radio 
emission. 

We recently found that at the center of the clusters A1835, 
A2029, and Ophiuchus, the dominant radio galaxy is surrounded 
by a diffuse low surface brightness mini-halo. We have studied 
the morphological and physical properties (i.e. length-scale, cen- 
tral brightness, emissivity) of these mini-halos by fitting their az- 
imuthally averaged radio brightness profile with an exponential. 
The method proposed here to derive the length-scale of halos 
and mini-halos seems promising because it is relatively indepen- 
dent of the sensitivity of the radio observations. The exponential 
model is attractive in its simplicity and involves a minimal set of 
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free parameters. While it cannot account for any local deviations 
from circular symmetry of the diffuse emission, this empirical 
method does provide a rough estimate of brightness and size of 
these sources. We compare the surface brightness profiles of the 
new mini-halos discovered in Paper I with data already available 
in the literature, both for mini-halos and halos. We find that the 
radio halos can have quite different length-scales but their emis- 
sivity is remarkably similar from one halo to another. This result 
could have important impUcations for theories of the origin of 
radio halos in clusters of galaxies but awaits confirmation from 
future, more sensitive, observations. In fact, if many faint ha- 
los have been missed by the current surveys, the average radio 



emissivity we found at 1.4 GHz, (J) ^ 10"*^ erg s"^cm"^Hz"\ 
should be considered an upper limit. In contrast, mini-halos span 
a wide range of radio emissivity. Some of them, like the Perseus 
mini-halos, are characterized by a radio emissivity which is more 
than 100 times greater than that of radio halos. On the other 
hand, the new mini-halos in cooling core clusters analyzed here, 
namely A2029, Ophiuchus, and A1835, have a radio emissivity 
which is much more typical of halos in merging clusters rather 
than to the mini-halos previously known. 

We discussed the exponential fit in comparison with the cur- 
rent theoretical models for the distributions of synchrotron elec- 
trons and magnetic fields in cluster of galaxies. We find that 
the exponential profile is very close to the expectations of these 
models if the magnetic field energy density roughly scales as 
the thermal gas density. In the framework of the equipartition 
model, the physical meaning of the Sr^. length-scale obtained by 
the exponential fit is that it marks the point at which the magnetic 
field strength is decreased to about half the value at the cluster 
center. Mini-halos would appear smaller because the magnetic 
field falls off more rapidly with radius in cooling core clusters, 
in agreement with the suggestions of Burns et al. (1992). 

Can we explain the larger dispersion in the emissivity of 
the mini-halos? We have very bright mini-halos, hke Perseus, 
which host a particularly bright compact radio source at their 
center. For these, we cannot rule out that part of the mini-halo 
radio emission is related to properties of the local intergalactic 



medium and part is correlated with the AGN activity of the cen- 
tral brightest galaxy, as suggested by the faint correlation, dis- 
cussed in Paper I, between the mini-halo and cD radio power. 
Indeed, the higher emissivity found in some mini-halos could be 
due to an extra-amount of energy supply from the strong AGN. 

On the other hand, we also observe low-surface brightness 
mini-halos, hke Ophiuchus, which host a faint AGN at their cen- 
ter. These mini-halos appear to be scaled down versions of the 
larger halos and they could be powered by similar merging pro- 
cesses. We note that, although cooling core clusters are generally 
considered relaxed system, when analyzed in detail they some- 
times reveal pecuUar X-ray features in the cluster center which 
may indicate a possible link between the mini-halo emission and 
some minor merger activity. Indeed, Burns et al. (2008) simu- 
lated the formation of both cool core and non-cool core clusters 
in the same numerical volume. These simulations confirmed that 
non-cool clusters are formed via major mergers early in their his- 
tory which destroyed the cool cores and left significant residual 
kinetic energy in the gas which might be used to power the ra- 
dio halos. In contrast, cool core clusters do not suffer any major 
mergers, thus preserving the central cool regions. However, the 
cool core clusters do experience regular smaller mergers which 
still inject energy in the intra-cluster medium, but more mod- 
estly than in the non-cool core clusters. These minor mergers 
might power the mini-halos. 
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Table 3. Radio information of halos and mini halos taken from the literature and reanalyzed in this work with the exponential model 
fit. 
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Col. 2: Type of diffuse emission contained (H=halo, MH=mini halo); 
Col. 6: FWHM of the circular Gaussian beam; 

Col. 10; Average radio emissivity over the volume of a sphere of radius 3rj, fc-corrected with a = 1. 
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